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Abstract 

Using a Corner Transfer Matrix approach, we compute the bipartite entan¬ 
glement Renyi entropy in the off-critical perturbations of non-unitary conformal 
minimal models realised by lattice spin chains Hamiltonians related to the Forrester 
Baxter RSOS models [30] in regime III. This allows to show on a set of explicit 
examples that the Renyi entropies for non-unitary theories rescale near criticality 
as the logarithm of the correlation length with a coefficient proportional to the 
effective central charge. This complements a similar result, recently established 
for the size rescaling at the critical point [23], showing the expected agreement 
of the two behaviours. We also compute the first subleading unusual correction 
to the scaling behaviour, showing that it is expressible in terms of expansions 
of various fractional powers of the correlation length, related to the differences 
A — Amin between the conformal dimensions of fields in the theory and the 
minimal conformal dimension. Finally, a few observations on the limit leading to 
the off-critical logarithmic minimal models of Pearce and Seaton [49] are put forward. 
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Introduction 


Entanglement is a very specific and intriguing feature of quantum systems [I]. It de¬ 
scribes truly quantum correlations between parts of a system. Besides the questions of 
principal nature that it arises in the interpretation and behaviour of quantum mechanics 
mm, it finds very interesting and promising applications in quantum information theory 
and quantum computing [3] , in condensed matter physics lai, as well as in the physics 
of black holes [7] and has even risen some interest in biological systems [8] . 

Many proposals have been formulated to quantify it (see for example [9] and refer¬ 
ences therein). In this paper we focus on bipartite systems, composed of a subsystem 
A and its complement A. In such case, the most convenient measure of entanglement is 
the Von Neumann Entropy restricted to A, the so-called Entanglement Entropy [1]. It 
is used as a measure of entanglement for pure quantum states. The behaviour of this 
quantity has been deeply studied in a wide range of systems from analytical, numeri¬ 
cal and experimental points of view. For example, interesting experimental protocols 
in many body systems have been recently proposed mm- It has been evaluated in 
different dimensions and in different regimes, but it is in 1 -|- 1 dimensions that it shows 
its most amazing mathematical properties. In particular, in the case of Integrable Mod¬ 
els, it reflects many mathematical features that otherwise would be difficult to probe. 
For example, the study of Entanglement Entropy in critical one-dimensional quantum 
systems is one of the most powerful techniques for the identification of the central charge 
c of the conformal field theory describing the low energy excitations of the model. This 
evaluation can be performed both analytically and numerically (e.g. DMRG |12] ) and 
usually requires less information than the study of the scaling of the ground state energy 
for the identification of the universality class. 

In the last decade there has been a particular focus on the evaluation of Entanglement 
Entropy in Integrable models, in particular in Conformal Field Theory (GET) |13l I14| . 
in their massive perturbations mm and in lattice spin chains [niiiiKisiisniiiiiiss]. 
Recently also non unitary models have been taken into account, both in the critical [23] 
and in the massive [ 23 | regimes. 

At first sight, non unitary theories might be considered just as non physical math¬ 
ematical curiosities, but there are many very interesting examples where they actually 
do play a physical role. For example, it is known that a strongly interacting 2D electron 
gas in a magnetic field produces edge modes described by GET minimal models, in what 
is known as Fractional Quantum Hall Effect (FQHE). In some cases it has been shown 
that these edge modes are described by non unitary GET |25) . In |26j . for example, the 
minimal model Ads ,5 has been considered. The non-unitarity arises from the fact that 
in this particular case the bulk is gapless (while in the ordinary case it is gapped) and 
then the edge can dissipate into the bulk. 

These excitations can be represented using generalised spin chains: the critical Fi¬ 
bonacci Chain (the “golden chain” |2^, minimal model ^^ 4 , 5 ), which represents the 
FQHE with filling u = ^ and its non-unitary generalisation to the Lee-Yang Model 
Ad 2,5 |2E] are among the most known chain representations. These critical chains all 
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belong to the wide class oilAq{sl{2)) invariant spin chains |29j in the universality class of 
conformal minimal models A4m,m', unitary or not. They can be further generalised in 
order to include their massive <hi 3 perturbations. The off-critical Hamiltonians of such 
chains have been recently obtained [30] and they are related to the RSOS models mm, 
through the usual connection between the time evolution operators of ID quantum spin 
chains and the row-to-row transfer matrix of 2D classical lattice models. 

The Hamiltonians of these chains are generically of the pseudo-hermitian type, i.e. 
they are not hermitian but have real eignevalues. In recent years they attracted a lot of 
interest in the description of many physical phenomena, ranging from optical effects to 
non-equilibrium systems, etc... Similarly, in Quantum Field Theory it has been widely 
discussed [33] how the scattering matrices of apparently non-hermitian theories can be 
physically well defined. In particular, using the definitions of [33], non-unitary theories 
can be classified as real or non real; in both cases the inner product is indefinite, but in 
real theories eigenenergies are real and the eigenvalues of the S matrix are pure phases, 
while in non real models the energies are not real and the S matrix has eigenvalues which 
are not pure phases. Notice that if the S matrix is a pure phase, there is a conserva¬ 
tion of probability through the scattering process; this conservation is a fundamental 
requirement for having a real, well-defined quantum physical theory. 

Our goal in this paper is to compute the bipartite Entanglement Entropy of the 
ground state for spin chain hamiltonians of the type introduced in m on an inhnite 
one dimesional lattice, taking as A the negative semi-axis and as A the positive one. 
Eor non-Hermitian hamiltonians this may sound ambiguous, as non Hermitian operators 
share the same right and left eigenvalues but have different right and left eigenvectors. 
However, one can argue [23] using PT symmetry and chiral factorization of GET that, at 
least at criticality, the left and right ground states do coincide, thus leading to a correct 
positive definition of the Entanglement Entropy for this state. 

This bipartite Entanglement Entropy calculation can be achieved by considering 
the Corner Transfer Matrix (CTM) introduced by Rodney Baxter [33] in integrable 
statistical systems on a 2D square lattice. This tool, originally used for the evaluation 
of partition functions and one point functions of the 2D models, has been extended to 
the evaluation of the Entanglement Entropy of the corresponding ID chains, following 
an approach developed in [Hi¬ 
lt is known that these chains, as well as their classical 2D counterparts, the class of 
RSOS^^s models, show different regimes with different physical behaviour. In the unitary 
case s = 1, Andrews, Baxter and Eorrester (ABF) [31] classified four regimes. There are 
two second order phase transitions: one between regimes I and H and the other between 
regimes HI and IV. The identification of universality classes at the transition points 
has been studied by Huse [33]. Approaching the HI-IV transition the RSOS^q models, 
for different r, are described by the universality class of unitary CFT minimal models 
In the I-H phase transition, instead, the approach to criticality is governed by 
the universality class of Z/, parafermionic CFTs [33] with L = r — 2. 

In the non unitary cases s 7 ^ 1 explored in |32) . the classification of regimes is more 
complicated. However, still there are regimes HI and IV (or X if r odd and s = 
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with a phase transition between them identified 1371138] with the universality class of 
non-unitary minimal models A4r-s,r- The identification of the universality class of the 
regime I-II phase transition is still an open problem in the s ^ 1 case. We shall not 
address this problem here and in the following we restrict our calculations to regime III. 

In the case of unitary ABF RSOS models, Franchini and De Luca [39] have recently 
used the CTM technique to compute the Entanglement Entropy, checking that in the 
phase transition III-IV it correctly gives the results expected in the conformal unitary 
minimal models Air-i,r while in the regime I-II phase transition, it gives those expected 
for the Zr _2 parafermionic CETs. 

In this paper we extend this Entanglement Entropy calculation to the regime III of 
s 7 ^ 1 RSOS models. The main goal of this investigation is to test on a concrete example 
of non unitary integrable theories where the calculations can be carried out exactly, 
the conjecture that the Entanglement Entropy scales logarithmically in the correlation 
length in a manner that parallels the scaling on a finite interval at criticality (for details 
of what we mean here see subsect. [^. 

The paper is organised as follows: 

• In section we recall a few notions on minimal GET models (unitary or not) and 
their characters, emphasising their relation with modular forms. 

• Then in section we introduce the basic concept and formulae of Renyi and Von 
Neumann Entanglement Entropies. We discuss in particular their behaviour for 
finite size of the subsystem at criticality and for infinte size but off-criticality. We 
consider in both cases the dominating logarithmic term, but also the form of the 
expected corrections. 

• In sectionj^we describe the Eorrester-Baxter (EB) RSOS models whose continuum 

scaling limit is described by the perturbed minimal model (where 

t > 0 is a perturbing parameter and the least relevant operator). 

• We briefly summarise in section [^ the Corner Transfer Matrix construction and re¬ 
view its use as a tool for the evaluation of partition functions and Renyi Entropy. In 
particular, specialising to the RSOS models, we emphasise the connection between 
the blocks in the partition functions on multi-sheeted surfaces with the conformal 
characters of the corresponding GET at criticality. 

• Applying the CTM tool to the EB RSOS model, in section we evaluate the 
Renyi Entropy for these theories. Our main results are the confirmation of the 
presence of the effective central charge in the leading scaling of the Entropy [23| . 
The difference from the unitary case where the usual central charge appears is 
due to the fact that the physical ground state and the conformal vacuum do not 
coincide in non unitary models. 

• In section we evaluate explicitly the most relevant corrections - nowadays tradi¬ 
tionally called unusual corrections following [40] - to the dominating logarithmic 
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behaviour. They exhibit power law decays with exponents given by a sort of effec¬ 
tive dimensions, i.e. the difference between the conformal dimensions of some of 
the relevant helds of the CFT universality class and the lowest (negative) conformal 
dimension that always appears in non unitary CFTs. 

• In section we briefly comment on the evaluation of Renyi Entropy in off-critical 
Logarithmic Minimal Models, which can be seen as particular limit realisations of 
the FB RSOS model. The observed absence of the double log behaviour is related 
to the destruction of the Jordan block structures responsible of the logarithmic 
behaviour as soon as the log-CFT are perturbed off-criticality. 

• Finally, in section we trace our conclusions and perspectives for future work. 


1 CFT minimal models and their characters 


We recall here some basic facts about conformal minimal models that we need in the 
following. The minimal models Mrn,m' are conformal theories whose Hilbert space is 
built up of two chiral parts each one composed of a hnite number of irreducible highest- 
weight representations (HWR) of Virasoro algebra at a given value of its central charge 
c. They are labelled by two coprime integers m, m' such that 2 < m < m! and have 
central charge 

6 (m — 

c = 1- - - 

mm' 


The (left) conformal families [‘ha,a'] are labelled by two integers a, a' running on the 
domain J = {[a, a') : 1 < a < m — 1, 1 < a' < m' — 1}. The conformal dimensions of 
the primary states |a, a^) of such families are given by 


^a,a' — 


(am' — a'my — (m — m') 
4mm' 


A2 


( 1 . 2 ) 


The Z 2 symmetry q/ = ^m-a,m'-a' is present in the whole series of models. 

The models are unitary for m' = m non unitary otherwise. In the non unitary 
case the state of lowest conformal dimension is not the conformal vacuum |0) = |4>i^i) 
but a different state |min) with negative conformal dimension (that can be proven always 
to exist and be unique) 


Amin — 


1 — (m — m') 
4mm' 


M2 


Correspondingly, an effective central charge can be defined 


(1.3) 


Ceff — C 24Aniin — 1 


6 

mm' 


(1.4) 


Notice that while the central charge c can be negative in non unitary models, the effective 
one Ceff is always positive. 
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The characters of the Virasoro HWRs can be written as 


Xa,a'{q) 



^k(km7n'-\-a7n'—a'm) _ ^(km-\-a)(km^^a') 


(1.5) 


where q = and (q')oo = n^i(l ~ They form a unitary representation of the 
modular group [H] PSL{2,Z) whose generators are 

T : T —>• r + 1 

c 1 

T 


fxa,a'{q) = Xa,a'iqe^n=e^^<^^P-^^)xa,a'{q) 
Sxa,a'{q) = Xa,a'iq)= s^J^^,Xb,b'{q) 


(b,b')£j 


where q = e ^ . 5 is the modular matrix 


= 2\l ^ (—+“ ^sin si 

’ V mm V m / 


Introducing the elliptic function 

OO 

fcSZ n=l 

they can equivalently be written as 


m \ 

( 111 

—-ab sm 

( TT —a h 

m / 

\ rn 


q^z-^){l 


/ f p ( m'a—a' 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 


iq) 


—a' m-\-mm' 2mm' 


-q^^E[-q 


^m'a-\-a'm-\-mm' ^2mm‘ 




The function E{z, q) is related to the first Jacobi theta function 


•')} 

(1.9) 


'dl{u-,p) = ^( — 1 )” 2p("-+2) g(2n+l)iu _ 2 pi “ 2p^‘^ cos2u +p^”')(l —p^"') 

nSZ n=l 

( 1 . 10 ) 

by the relation 

E{e^^^,q) =iq-h^'^Mw,qh (TH) 

The first Jacobi theta function, under a modular S tranformation, behaves as 

di(z,e“^) = —iy/ire^^ 'di{TZ]e™'^) (1-12) 

These formulae will be useful later in sect, [^for the computation of the Renyi Entropy. 

The minimal models can be perturbed off-criticality by picking up combinations 
of their relevant fields, resulting in super-renormalizable theories. In particular, the 
perturbation by the least relevant field <bi ^3 preserves the integrability of the model off- 
criticality. It is this integrable perturbation that represents the scaling limit of the RSOS 
models in the off-critical vicinity of the phase transition between regimes III and IV. 
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2 Renyi Entanglement Entropies 


Consider a bipartite quantum system whose Hilbert space is = 1-La ® Ha- If 
system is in a pure state |'I') (i.e. one with density matrix p = |'I')('b|), the entanglement 
between the two parts A and its complement A can be described by the family of Renyi 
Entanglement Entropies Sn {n G R>o) 

S^A = V— log ^^aPa (2-1) 

where pA = 1>^aP is the reduced density matrix of the subsystem A. The Von Neumann 
Entanglement Entropy S = lim = —tr^ ipA^ogpA) is a special case of the Renyi 

n— 

Entropies. 

In |13[ fH] it has been shown that the Renyi Entropy for a bipartite one-dimensional 
quantum system scales as 

(n) C n -|- 1 Z (r, r,\ 

Sa ^ a -^og- (2.2) 

on e 

where I is the size of the subsystem A, and e is some ultraviolet cut-off. This scaling holds 
when the size of the whole system is much larger than the size I of subsystem A and when 
the correlation length ^ is much larger than 1. This condition on the correlation length is 
equivalent to assume that the system is approaching criticality with a universality class 
identified by a GET with central charge c. In |40j the subleading contributions to this 
behaviour were estimated and look like 

sin) ^ l_ ^ ^ ^in)i-2xln ^ y{n)i-2x ^ _ _ (2.3) 

on e 

The constants and are non-universal, but the exponent x is, and has 

to be identified with the scale dimension x = A -|- A of some operator <I> of the GET 
universality class having the lowest conformal dimension among those concurring to the 
corrections. Higher order corrections are expected and they are due to other primary 
fields and to descendents. For n > 1 the term is subleading with respect to the 
l- 2 x/n Viceversa for n < 1. Taking into account both terms is necessary to ensure 
a smooth limit for n —?• 1. 

On the other hand, if the correlation length ^ is finite and the size Z —)• oo we have 
|42| a similar expansion 

~ + ... (2.4) 

which represents the entanglement between two infinite parts of an infinite system. In 
this case the size I of the subsystem A has been replaced by the correlation length. This 
is due to the fact that the leading term of the Entropy is given by the smallest physical 
length which plays a role in the system. The factor of 2 difference in the leading terms 
and in the corrections is simply explained by the fact that the number of boundary 
points dividing A from ^ is 2 in the first case and 1 in the second. The non-universal 
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coefficients of the c orrec tions in general can be different from the 

a(”'),&(”■), of eq.(2.3). Also, the universal exponent h can differ from the x above 

[2ll[39]. 

Recently [2H] this evaluation has been extended in order to take into account non 


unitary systems. The Hamiltonian operator H describing a CFT system of Hilbert space 
7i put on a cylinder of circumference R is given by 


H = — (Lo + Lo- — 

12 


(2.5) 


where Ln,Ln are Virasoro algebra generators and c is their central charge. In unitary 
CFT the conformal vacuum |0) (defined as the unique state such that Ln\0) = 0, Vn > 
— 1) is the physical ground state, i.e. V|V’) G R '■ (0|i7|0) < {^p\H\^p)). Thus the ground 
state energy will scale as Eq = —^ |l3l 01] . 

In non unitary CFT instead, there is at least one primary state |min) with negative 
conformal weight Amin ((min|Lo|iiiin) = Amin < 0). Obviously, this state and not 
the conformal vacuum | 0 ) has the lowest possible energy on the cylinder and is the 
true ground state of the theory. For this reason the ground state energy scales as 
F'min = — , with the so called effective central charge Cefr = c — 24Amin > 0 instead 

of the central charge c [45| . The difference between the conformal vacuum and the 
physical ground state plays an important role in the definition and calculation of the 
Entanglement Entropy. For a detailed discussion, see [23] . As a result, in non-unitary 
CFT the Renyi Entropy scales as 

Q{n) Ceff n -|- 1 / ("O 

^ — -log - ( 2 . 6 ) 

bn € 

Although this result seems natural, its proof in a generic non unitary CET is far from 
trivial and can be acheived by a careful analysis of how to correctly define twist operators 
in such case through an orbifold approach [23] . 

In the case of infinite subsystem, one is tempted to replace ^ by the correlation 
length, as in the unitary cases. In this paper our aim is to verify that this conjecture is 
correct in a concrete lattice realization of an off-critical non-unitary model. We use the 
RSOS non-unitary models of Forrester Baxter in regime HI as lattice regularizations of 
the perturbed CFT A4m,m' + (where t > 0 is a perturbing parameter). 

In other words, we try in this paper to reproduce a formula of the type 

5^^ ~ (2.7) 

calculating it from the exact expression for the Renyi Entanglement Entropy obtainable 
by a Corner Transfer Matrix approach combined with the evaluation of the correlation 
length ^ on the lattice. We expect that the exponent h is affected by the non unitarity 
in a way similar to the change from c to Cefr in the leading logarithmic term. 
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Figure 1: Square tile 


3 RSOS FB Models 


The Restricted-Solid-On-Solid RSOSr,s model [SU [32] is defined on a square lattice, for 
each pair of coprime positive integers r, s such that r > 3 and 1 < s < r — 1. On each site 
i £ 7? there is a variable, called local height, li that takes values = — 1. Local 

heights of two neighbouring sites % and g are restricted to differ by one — = 1. Local 

heights can be thought as encoded on a A^-i Dynkin diagram. Each node represents 
a possible value that £i can take and it is linked to the possible values at neighbouring 
sites. This allows to generalise the model, following Pasquier |46j . to other simply laced 
A, D, E Dynkin diagrams. For the scope of the present paper, however, we focus on the 
A case only. The RSOS models belong to the wide class of 2D classical lattice model 
of IRF (Interaction Round a Face) type. In IRF models, the interaction is on nearest 
neighbour, so that one can define local Boltzmann weights that depend on the four sites 
i,j, k, I around a tile (see fig. 




(3.1) 


where /3 is the inverse temperature, is the energy of the configuration of the four 
vertices. 

For the RSOS models, the Boltzmann weights have been calculated from Yang- 
Baxter equation, that must be satisfied by any integrable lattice model, in mm- 
They depend on a spectral parameter u, measuring the anisotropy of the lattice, on the 
crossing parameter 

A = -TT (3.2) 

r 

ruling their behaviour when the lattice is rotated by 


W 


ii ij 


s(^^A)s(4A) 

si£jX)s{£iX) 


£k 


;! X-u\ 


(3.3) 


and on a temperature-like parameter p, (—1 < p < 1), such that t = p^ is measuring the 
departure from criticality. The system is critical for p = 0. The Boltzmann weights turn 
out to be expressible in terms of elliptic theta functions for which the parameter p plays 
the role of the nome. Here and below s{u) = '&i{u-,p), the first Jacobi theta function 
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(1.10). The non-zero Boltzmann weights can be put in the following form (in the so 
called symmetric gauge) 


w 

£±1 

£ 

£ 

w 


£±l 

£ 

w ( 

' £ 
£±l 

£±l 

£ 


u = 


u = 


s(A — u) 
s{X) 

^/s{{£ =F l)A)s((£ ± 1)A) s{u) 


s{iX) 


s{X) 


(3.4) 


u = 


s(£X ± u) 
s(£X) 


The periodicity and modular properties of the Boltzmann weights determine the funda¬ 
mental ranges of the parameters. The classification of possible regimes of the models is 
quite complicated |32j but in this paper we concentrate on regime III (0 < p < 1 and 
0 < u < A) only, leaving the others for future investigation. 

For s = 1 the model, often denoted simply RSOS^, has been investigated by Andrews, 
Baxter and Forrester in m- The Renyi Entanglement Entropy for the cooresponding 
spin chains has been computed in |39j , where it is shown that in the regime III approach¬ 


ing the critical point p = 0 it scales like (2.4) with the central charge of unitary minimal 
models In the regime II transition, instead, it scales as the parafermions 

Here we are interested in the generalizations of these results for other values of s. 
The spin chains Hamiltonians corresponding to the RSOSr,s models have been recently 
written in m- As the notation for them is quite complicated and not relevant for the 
following, we do not write them here and invite the interested reader to refer to the 
paper m for an exhaustive presentation. 


4 The Corner Transfer Matrix approach 


Eor spin chains on an inhnite ID lattice, we consider the Renyi Entanglement Entropy 
between two semi-infinite halves, the negative one conventionally called A while the 
positive is A. For chains related to a 2D classical lattice model of IRF (Interaction 
Round a Face) type, an efficient method to compute the reduced density matrix pA, as 
proposed in is to make use of the Corner Transfer Matrix approach [33]. The FB 
RSOS models are of this kind, therefore we adopt this approach here. 

Consider the 2D “diamond shaped” lattice of hgj^ and divide it in four quadrants. 
The central site is denoted by 0. In the lower-left quadrant A introduce an operator (see 

fig- ^ 



{u) 


EH”' 

DeA 


£i 

£^ 



(4.1) 


where £ = {£o,£i,£2, ■■■,£n) and £' = {£'q, £'i, £' 2 -, ■■■,£' j^) are vectors collecting all the 
variables along the two inner boundaries, i.e. on the sites denoted by o in fig. The 
sum is performed over all possible values of £i on internal sites (signed by a black dot • in 
fig. [§. The variables at the outer boundary sites are assigned fixed values determining 
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Figure 2: The 4 corner transfer matrix operators 


a boundary condition uniquely. The product is over all tiles □ of the quadrant. Notice 
the Iq = fg obvious constraint on the central site 0. 

Analogously, define in the other quadrants the operators and 

Organizing the product of VF’s to be performed diagonally (thus the diamond shape of 
the lattice), a corner transfer matrix like may be expressed in terms of smaller 

corner transfer matrices and A^^, (dehned for reduced quadrants). This 

recursion relation allows, in principle, the iterative calculation of the corner transfer 
matrix for any lattice quadrant of finite size. 

The partition function can be expressed in terms of the four CTM operators as 

= tr(A(^)B(^)c(^)D(^)) (4.2) 

There are other interesting quantities that can be computed form the CTM method, like 
e.g. the height probabilities. Here we are interested in the calculation of the reduced 
density matrix pA relative to a dominion A coinciding with the negative axis of the 
horiziontal direction of the infinite lattice that results taking the thermodynamic limit 
A —)■ oo of the present construction. Following the unnormalised reduced density 
matrix on the subsystem A can be written a^ 

Qa,W = (ABCD)££/ (4.3) 

^Please notice the difference of the symbols q (unnormalised density matrix) and p (normalised density 
matrix) 


11 



























































where A = limjv->-oo and analogously for B, C, D. Obviously tr aqa = -2^1 is the 
partition function in the thermodynamic limit 


= lim 

N^OO 


(4.4) 


To normalise, one has to divide by the partition function 

_ (ABCD),,, 

PA,ee' — 


(4.5) 


For the Renyi Entanglement Entropy we need to compute ti ap'a- Defining the higher 
genus partition functions iiAP^ = -2^n one sees that 


truPu = 


Zn 

Z’^ 


and the expression for the Renyi Entropy is 



1 

1 — n 


log 




(4.6) 


(4.7) 


As shown by Baxter [33], the Yang-Baxter equation that has to be satisfied by 
the Boltzamnn weights of any integrable model, implies that the four CTM operators, 
for A —)■ oo, commute each other for different values of their parameters and their 
eigenvalues, up to a common factor, accommodate into a diagonal matrix of the general 
form 

4"/,, = (4.8) 

The diagonalisation can be performed along the lines illustrated in |34] . In particular, 
for FB RSOS models it has been performed in the original paper [32|. It is simplified if 
one introduces new variables x, y that in regime III are defined as 


y 


47r^ 

glogp 


47r^ s s 


(4.9) 


Define 


N 


m = Y /i;7(4,4+i,4+2) 


k=l 


wher^ 


7(^,£±1,4 = ^ 

In regime III we have 


4 

r 


7(^±l,f’,^Tl) = 2 


(4.10) 


(4.11) 


4>(4 = ^ I+ (4 - 4+i) 

k=l 



^ [x\ denotes the integer part of x. 


(4.12) 
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The second term of the summand is absent in the s = 1 case simply because is 

always 0. Consider now a system with size 2N + 3. The unormalised reduced density 
matrix in regime III, as computed in [32], is given by 




, = (ABCD)i^) = E (x^°,y) 


(4.13) 


where Iq is the central height and £ = {£q, ... ,£Ar+i) is a vector of all local heights from 


the central to the boundary £n+i- Here E{-, •) is the modular form defined in (1.8). 


5 Renyi Entropy in FB RSOS models 

As disccussed in section to compute the Renyi Entropy we need to consider not only 
QA but also its re-th powers q'\. 

Let us set the boundary condition to = 6 and £n+i = c = b ± 1. The n-th power 
of the reduced density matrix is given by 


The trace can be performed summing over all allowed configurations £ 

zr> = 

I 

r—1 

a=l 


(5.1) 


(5.2) 


where D]\f{a, b, c]q)= ^ with £o = a, £n = b and £]\f+i = c. The symbol 

means that the sum is restricted to even or odd only values of a accordingly to the 
parity of the boundary conditions. In particular notice the the central height a and the 
boundary height b have to be compatible with the requirement \£i — £i+i\ = 1. 

The limit A —)> oo can be performed, keeping track of the boundary condition ^ 
£n+i'- 

kN . ,1 6(6-1) {fe-l)6 

lim q 2 DM{a,b,b+l]q) = ——q 4 2 F{a,b-k-,q) 

N^oo {q)oo 

with k = if < £n+i or 


(5.3) 


M. , ,1 6(6+1) fc(6+l) 

lim g 2 E)jY(a, 6 + 1, 6; g) = —— q 4 2 E[a,b-k;q) (5.4) 

N^oo (g)oo 


with k= [f\ = if 

The function F is defined as 


a(a —1) 

F{a,d;q) = q~i~ 


ad 

2 


= Q 


E ^_qrd+('>'-a.){r-s) ^^2r{r-s)^ _ ^adj^ ^_^rd+(r-a){r+s) ^ ^r{r-s) 

(a-l) ad (g) 


q 


-^+A„ 


-Xda{q) 


(5.5) 


13 

















where the Virasoro characters XdaiQ) (see eq.(1.9)) are taken here for m = r — s and 
m! = r. Thus 


kN r /L J'\ q(q. —1) _od _|_ _c_A 

lim q+ 2 D]si{a]q)b,d=b-k = ' ’q ^ 2 +24 '^^Xda{q) 


N^OO 


where 


f+{b,d) = 
f-{b,d) = 


b{b-l) b{k-l) 

4 2 

6(6 + 1) fc(6 + l) 


(5.6) 


(5.7) 


do not depend on a. 

Since the physical (normalised) density operator is given by 


PA = 


QA 


ABCD 


TiaQa Tr(ABCD) 
the quantity in which we are interested is 


i n 


where 


r —1 


Zn= E{x'^,yy 




X 


Xda (a;^”) 


a=l 


(5.8) 


(5.9) 


(5.10) 


since the factor (/i^ + appears in both the numerator and the denominator and 
cancels out. 

Notice that the critical point is reached for p —)• 0, or x —)■ 1. In order to catch the 


critical behaviour we transform (5.10) into a new expression which is more suitable for 


expansion near p = 0. First of all we can use the relations (1.11) and (1.12) for the 
function E 


r—1 

/ TTaS /—\^ { 2r?\ 

Wn = ^ [ — ^Vp) Xda ) 

a=l 


(5.11) 


where, again, Wn is proportional to Zn and a common factor with Z^ has been dropped 
so that 

, n kVn , . 


Furthermore we can perform a modular S transformation (1.6) on the conformal 
character 

-b,b '. 


Xa,a'{Q)= '^aA'Xb^b'iq) 

(b,b')ej 


(5.13) 
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^ s 27r2 r_ i ^ ^ P ^ 

Since = e V logpry^ niodular 5-transform is cc = (x^) = e 2 s = 

(recall that, in regime III, p > 0). Similarly, (x^"') = {p^)^ = We have 


r—1 


a=l d'a 


w„ = Y.' E'’> (^. Vp)"st'x,iv (-i 


(5.14) 


which is suitable for a p —)■ 0 expansion. 

Denoting f) = (d', o') G ^7 as an index which spans the Kac table we have 


y^n = ^xtj fli(n,p) 


l)ej 


— Xmin ( ^ 


^ 1 /min p) 1 + 


f) 7 ^min Xmin I ^ 


XI, fi,(n,p) 

k] /min (n,p) 


(5.15) 


where 


r—1 


a=l 


„ / s V—T* / TTdS \ ^rl'n' 

h {n,P) [ — ,Vp) 


(5.16) 


and 1) = min refers to the primary field with the lowest conformal dimension Amin. 
Taking the logarithm, we have 


1 1 

logTTT^ = log 


Xmin I W " 




+ log 


/min (n, p) 


Xmin {ojT /min(l,p)” 

/„(n,p) 


+ log 1+ X] 


/min (n,p) 


— n 


log 1+ X] 


ll/min Xmin 

Xl,(w) /f,(l,p) 


[j^min (l-iP) 


(5.17) 


Expanding Xmin and /min of the first two terms of the equation above near p, cu = 0 we 
obtain the leading scaling and the constant coefficient for the Renyi Entropy 


q(«) _ 


+ 


Ceff n + 1 

24 n 


logo; Til’ll 


h{n,p) 

tj^min Xmin (w") 


rvi;>°eh+ E 

y f)7i:mi 

-login- W 


Xmin (^) /min (iiP) 


(5.18) 

(5.19) 
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The constant is given by 


r—1 



which is well defined in the n —>■ 1 limit. 


(5.20) 


6 “Unusual” corrections 


The next step is the evaluation of power law corrections to the logarithmic scaling of the 
Renyi Entropy. First, notice that Xh ^ Xmin for cj —)• 0. In this regime the argument of 
the logarithms in the second and third line of (5.18) is close to 1 and then it can be Taylor 
expanded. Taking into account only the most relevant contribution and assuming, to fix 
ideas, that n > 1, we have 


c,(n) _ _ 

“ 24 


Ceff n + 1 


n 


logw + i^") + 


E 


with 




E - n nas c’t 

sm 

^ a 

l-nJ2 sin” 


( 6 . 1 ) 


( 6 . 2 ) 


and the most relevant contribution is given by the second smallest conformal dimension 
among those appearing in the expansion of |5.18 denoted Ai in the following. Taking 
into account only this correction we have 


Q{n) Ceff n + 1 aM , I 

—-tTa - loguj + A'' ’ + B] 'uj -\ - 

^ 24 n ^ 


In the case n < 1 one can proceed similarly, but now the first correction will be 

sin) ^ 1 ^ J^(n) ^ ^ . 

^ 24 n ^ 


with 


fb'(n) _ 


n 


Esin- 


ras on 
r ^da 


1 - n E sin ^5, 


= nB. 

min 0 

da 


(1) 


(6.3) 


(6.4) 


(6.5) 


In order to interpret these results from a physical point of view, i.e. to have an expression 
for the entropy in terms of the correlation length (or in term of the mass), we need a 
relation between the parameter ui = and the correlation length In other words, 
we need to know the critical exponent ly in regime III: 


m 


= 


( 6 . 6 ) 
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Using perturbative CFT it is possible to evaluate the critical exponent v. In the 
continuum limit, the model is described by the perturbed CFT: 

S = S'cFT + t J (fx d>i, 3 (a:) (6.7) 

where |t| = measures the departure from criticality. A simple dimensional analysis 
of this action tells us that v = 2(i-Ai 3) ~ 37 minimal model A4r-s,r at the 

transition between regime III and IV. 

This value for the critical exponent can also be extracted extending previous results 
m to the Forrester Baxter models. Taking into account the necessary modifications, 
the calculation can be carried out along the same lines of 133 to get 

e”? = k'{p^) (6.8) 


where k'{q) is the conjugate elliptic modulus for the elliptic nome q 


t'(<i)=n(i 


i=l 


1 — q 
+ 9 


21-1 


Expanding (6.8) around p = 0 we get 

r' = f f + o(p9‘^) 

357 

in perfect agreement with the perturbative CFT prediction. 

Using (4.9) and (6.10) we have u = (8^)“^ + • • • which gives 

^ log^ + ^ 

where constants and are just trivial rescaling of and 

^(n) ^ i(n) + ^!i±llog2 

4 n 


(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 


From the equation (6.12), it is immediately possible to read the correction to the loga- 

- 2A 

rithmic scaling. For unitary theories the correction is expected to scale as ^ « where A 


is the conformal dimension of the field related to the correction [39] . In the non unitary 
case this term is affected by the fact that the ground state is no more the conformal vac¬ 
uum |0) but another state |min) = <hmin(0,0)|0). The formula Ai — Amin shows that the 
presence of a non trivial ground state modifies the scaling of the correction. It is known 
|23j that this feature affects also the logarithmic scaling of the entropy in non-unitary 
models, where the central charge c is replaced by the effective one Ceg. Similarly, it is 
not surprising that the conformal dimensions are replaced by some sort of “effective” 
dimensions A — Amin- Notice that the unitary case can be immediately recovered in 
(6.12) by setting Amin = 0. 
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7 A comment on off-critical Logarithmic Minimal Models 


An interesting feature of Forrester Baxter RSOS models is that, for a particular choice 
of parameters r and s, they provide a lattice realisation of Logarithmic Minimal models 
j48| and their off-critical thermal perturbations |49j . In particular it has been shown 
that in the limit r, s —)> oo, keeping the ratio r/s fixed to some rational number R/S, 
the underlying model becomes the so called logarithmic minimal model CAiR-s,R and 
its off-critical thermal perturbation CM.r-s,r + t^i ,3 ■ 

Taking such a limit for the entropy, the result is not affected from the functional 
point of view: it maintains the same structure and the effective central charge 


CpflF = lim 

r,s—>oo 


1 - 


6 


r(r — s) 


(7.1) 


is identically 1 for any choice of the ratio R/S 


= 


1 n -b 1 
12 n 


logC + -b 


(7.2) 


Here = lim . When taking the limit r, s 

r,s^oo 


oo the S modular matrix becomes 


ill defined because the prefactor 


tends to zero in this limit. While 


(n) 

this feature does not affect the limit of the coefficent ' (6.2) - as 5 appears with 
the same power both at the numerator and at the denominator - it implies some extra 
attention for the limit of the constant (5.20). For this reason, we need to modify 
the definition of the Renyi Entropy for Logarithmic models multiplying the partition 
function by r in a sort of renormalisation procedure 


^(n) 


lim log-) 

r,5^oo y T J 


(7.3) 


This multiplication can be seen in the same spirit of defyning generalised order param¬ 
eters in [39], while a better understanding of this renormalistion of the Entropy is still 
missing and goes far beyond the aim of this work. In any case, it only affects the nu¬ 
merical value of the non-universal constant and not the functional shape of the ^ 
dependence of the entropy. 

This ^ dependence result (7.2) disagrees with the general prediction for the Renyi 
Entropy I dependence in logarithmic CFT |23j , where a double logarithmic (log log) term 
is expected. The difference is due to the fact that the logarithmic feature of the system 
is related to the presence of non diagonalisable Jordan blocks in the Hamiltonian. It has 
been noticed in |39| that such non diagonalisability feature disappears as soon as the 
system is perturbed thermally out of the critical point. Thus the so called off-critical 
logarithmic minimal models are not really logarithmic and thus we should not expect 
double log term in the computation of the entropy. In the case of logarithmic minimal 
models, then, we expect that, while the dependence in the subsystem size I at criticality 
rescales with a double logarithmic term, this is absent in the rescaling off-criticality in 
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^ while approaching the logarithmic critical point. The two behaviours in this case are 
different, while in the usual minimal models (unitary or not) they show a similar (single 
log rescaling) pattern. 

8 Conclusions and outlook 

We have computed the Entanglement Renyi Entropy for Eorrester-Baxter RSOS models 
in the regime III, a set of lattice systems whose continuum limit is described by non¬ 
unitary Minimal Conformal Models perturbed by their $ 1^3 operator. Our evaluation 
focused on the scaling with the correlation length ^ of the Entropy for a infinite, bipartite, 
quantum system. The computation of Entanglement Entropy in non unitary models has 
been recently addressed in |23j and [23] where a logarithmic scaling with the size I of the 
subsystem has been found. In particular it has been shown in [23] that the coefficient 
of the logarithm is given by the effective central charge Ceg. Eor unitary theories the 
computation of the Entropy for a finite critical system [I < 00 , m~^ = (, = 00 ) can 
be easily translated [IT] in the infinite subsystem with a small but non zero mass gap 
(Z = 00 , m~^ = ^ < 00 ) using arguments similar to the Zamolodchikov c-theorem. 
In the non unitary case such a translation cannot be performed a priori, since some 
assumptions about the positivity of certain correlation functions are no more valid in 
the non unitary case. Eor this reason, the scaling we have found, where the proper length 
I is replaced by the correlation length ^ gives a strong hint about the renormalisation 
flux properties in the non unitary case. A proof of a sort of Cefj-theorem in the non 
unitary case goes far beyond the aim of this paper. Nevertheless, our result supports 
evidence, in a non trivial set of non unitary models, that the ^ scaling follows the same 
logarithmic law of the I scaling, exactly like in the unitary models. 

We have also studied the power law corrections to the logarithmic scale of the En¬ 
tropy. The interesting part is not in the coefficients of the expansion, that are non 
universal, but in the exponent of the power of In the unitary case the expansion or¬ 
ganises into many series expansions in terms of powers of or with f) labelling 

the various conformal primary fields. In the non unitary case we find that the series 
expansions are in terms of qj- ^(A(,-Anjin)/n other words, also the conformal 

dimensions here take an “effective” value shifted by Amin, exactly like the central charge 
is replaced by Ceg. 

We have also briefly considered the limiting case of the off-critical Logarithmic Min¬ 
imal Models. Even if a double logarithmic term has been expected from the literature 
[23], where the I rescaling has been considered, we did not find here such a term; the 
scaling in ^ behaves exactly as in the non-logarithmic case. This discrepancy is due 
to the fact that the perturbation outside the critical point destroys the Jordan blocks 
responsible for the creation of the logarithmic features, as already pointed out in |49| . 

It their seminal work on the EB RSOS model [32|, Eorrester and Baxter classified 
many regimes, while we have restricted our analysis to regime III only. While it is 
known [371 EH] that the regime III is a lattice realisation of the perturbed Minimal 
Models Air-s,r + ‘hi, 3 ) the physical interpretation of the other regimes has still to be 
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explored. In the unitary ABF RSOS model m, the number of regimes is smaller (just 
four, compared to ten in the FB case) and two kind of universality classes have been 
identified. The critical line between the Regime III and IV can be described by the 
Unitary Minimal Models A4r-i,r, while a unitary parafermionic CFT Zr _2 underlies the 
critical line between regimes I and II. It would be very interesting to try to understand if 
some critical line of the FB RSOS model can be classified using, maybe, the non unitary 
parafermionc CFT introduced in [5^ or if they obey instead some other classification 
scheme. The evaluation of Entanglement Entropy has been demonstrated to be a pow¬ 
erful tool for the identification of universality classes and it could be a valid help also in 
the study of the unknown regimes of the FB RSOS model. 

Of course, further generalizations of these results to the lattice realisations of higher 
coset conformal filed theories perturbed by their opertors, or even for models based 
on cosets of higher algebras, like Wn-algebra based series, including Pasquier like gener¬ 
alizations and their dilute versions, etc... can all be the subject of future investigation. 
As non unitary coset models as critical points of fusions of FB RSOS models are under 
construction right in these days m, the exploration of Renyi Entanglement Entropies 
for this set of models is viable. 

The bipartite Entanglement Entropy can give a lot of useful information, but a more 
complete understanding of entanglement needs also the knowledge of the dependence 
on the size of the system, as well as its behaviour in finite intervals. We know how 
to deal with these problems in GET m and, with a form factor expansion, also in 
off-critical integrable quantum field theories m- However, in integrable spin chains, a 
useable procedure to compute entanglement entropies on a finite interval is still lacking 
in general. This would be a serious progress in our understanding of entanglement in 
one-dimensional systems. 

Also, the interest in entanglement measurements with more than two subsystems 
introduces new quantites to be evaluated, like negativity [52l [53]. To find a way to 
compute these quantites in a consistent lattice integrable approach is a challenge for 
future research. 
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